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We report on the first axion search results from the Cryogenic Dark Matter Search (CDMS) 
experiment at the Soudan Underground Laboratory. An energy threshold of 2keV for electron- 
recoil events allows a search for possible solar axion conversion into photons or local Galactic axion 
conversion into electrons in the germanium crystal detectors. The solar axion search sets an upper 
limit on the Primakov coupling g all of 2.4xl0 -9 GeV -1 at the 95% confidence level for an axion 
mass less than O.lkeV/c 2 . This limit benefits from the first precise measurement of the absolute 
crystal plane orientations in this type of experiment. The Galactic axion search analysis sets a 
world-leading experimental upper limit on the axio-electric coupling g a ee of 1.4xl0 -12 at the 90% 
confidence level for an axion mass of 2.5keV/c 2 . This analysis excludes an interpretation of the 
DAMA annual modulation result in terms of Galactic axion interactions for axion masses above 
1.4keV/c 2 . 

PACS numbers: 14.80.Mz, 29.40.-n, 95.35.+d, 95.30.Cq, 95.30.-k, 85.25.Oj, 29.40.Wk 



The axion has been postulated to solve the strong 
CP problem in quantum chromodynamics. The break- 
ing of Peccei-Quinn UYl) symmetry leaves a pseudo- 
Goldstone boson field interpreted as the axion. Al- 
though the original Peccei-Quinn axion model has been 
ruled out [2|, l3j, "invisible" axion models allow a wide 
range of axion masses and axion-matter couplings 
Astrophysical observations are currently the best strat- 
egy to search for these invisible axions The interior 
of stars is expected to be a powerful source of axions due 
to the high abundance of photons and strong electromag- 
netic fields, which may convert photons into axions. The 
non-thermal axion production mechanism in the early 
universe provides a cold dark matter candidate. 

Here we report on the first axion search results from 
the Cryogenic Dark Matter Search (CDMS) experiment. 



The CDMS collaboration operates a total of 19 Ge 
(~250g each) and 11 Si (~100g each) crystal detectors 
at ~ 40 mK in the Soudan Underground Laboratory. The 
detectors are designed to read out both ionization and 
phonon signals from recoil events [6J . The ratio of ion- 
ization to phonon energy, the ionization yield, enables 
discrimination of nuclear from electron recoils. The de- 
tails of the detector structure and operation can be found 
in reference [2j. We report on data from two run peri- 
ods between October 2006 and July 2007, also used for 
the WIMP-search analysis described in [8[ . This analysis 
modifies our event selection criteria to focus on electron- 
recoil events, lowering the energy threshold to 2 keV for 
such events and giving a net Ge exposure of 443.2 kg-days 
before cuts. 

The flux of solar axions at the Earth can be estimated 
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assuming the standard solar model and a coupling 
to the C(keV) black body photons in the core region of 
the Sun. For axion masses <C 1 keV/c 2 , photon-axion 
conversion creates a flux of O(keV) axions at the Earth. 
The solar axion flux at the Earth is given by 13, II | : 
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where E a is the energy of the axion in keV and <? a77 is 
the axion-photon coupling constant. 

The axion-photon coupling to the nuclear Coulomb 
field in the detectors converts axions back into photons 
of the same energy (Primakov effect). Coherent Bragg 
diffraction produces a strong correlation between inci- 
dent beam direction and conversion probability, provid- 
ing a unique signature of solar axions. The expected 
event rate can be computed as a function of energy and 
the orientation of the crystal relative to the location of 
the Sun [lj|. It is given as a function of observed pho- 
ton energy E for a given axion momentum transfer q and 
scattering angle 9 [13J, by: 
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where W is a detector energy resolution function. The 
Fourier transform of the electric field in a crystal is 
given as F(q) = k 2 J d?x<j)(x)e n \ which depends on 

<t>{%) = J2t 47r\s-x t \ e ~ lJLllA = Y,G n Ge lGS , where k is 
the photon momentum, e is the elementary charge, Xi 
is the position of a germanium atom in the lattice, r is 
the screening length of the atomic electric field, Z=32 
for germanium, and G is a reciprocal lattice vector. The 
structure coefficients uq (defined in [13l ]) account for the 
face-centered-cubic structure of Ge. The Bragg condition 
(q = G) can be expressed in terms of the axion energy as 
E a = hc\G\ 2 1 (2u ■ G), where u is a unit vector directed 
towards the Sun. 

The expected event rate is calculated based on an ac- 
curate measurement of the orientation of each detector 
with respect to the position of the Sun. We took the 
specific geometry of the experiment, the live-time dur- 
ing data taking and the seasonal modulation of the so- 
lar axion flux due to the changing distance between the 
Sun and the Earth into account. The geodesic loca- 
tion of the Soudan Underground Laboratory is latitude 
47.815°N, longitude 92.237°W and altitude 210 m below 
sea level. The geodesic north of the CDMS experimental 
cavern was measured in 1999 by the Fermilab Alignment 
Group 



14j . A line connecting two survey points along 



the central axis of the cavern was found to be 0.165°E 
from true north. By extension, the main horizontal axis 
of the CDMS cryostat was found to be 0.860 ± 0.018°E 
from true north. 
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FIG. 1: Time and energy dependence of the expected solar 
axion conversion rate in a Ge detector for g all = 10 -8 GeV -1 . 



Within the cryostat the 30 CDMS detectors are 
mounted in five towers of six detectors each. The ver- 
tical axis of each tower is aligned with the [001] axis of 
the detectors. The (110) axis that defines the major flat 
on each substrate is rotated with respect to its neighbors 
above and below, such that the detectors form a helix 
within each tower. The uncertainty in the absolute az- 
imuth orientation of the crystal planes is dominated by 
an estimated 3° uncertainty in the exact angular posi- 
tion of the tower axes with respect to the central axis of 
the cryostat. The uncertainty of the zenith angle mea- 
surement was estimated to be less than 1°. In Fig. Q] we 
present the predicted event rate in a germanium detector 
for an assumed coupling of g a ^y — 10" 8 GeV" 1 . 

In order to sample pure axion interaction candidate 
events, we require that an event (> 3cr above mean noise) 
was recorded in one and only one detector; i.e., a single 
scatter. To make sure the selected events are not due to 
residual cosmic ray interactions, they are required not to 
be coincident in time with activity in the veto shield sur- 
rounding the apparatus. Candidate events are selected 
within the ±2ct region of the electron-recoil distribution 
in ionization yield. Data sets taken within 3 days after 
neutron calibrations are not considered in order to avoid 
high gamma rates due to activation. The detection effi- 
ciency is dominated by rejection of events with an ioniza- 
tion signal in a detector annular guard electrode, covering 
15% of the detector volume. The acceptance of our anal- 
ysis cuts for single-scatter electron recoils was measured 
as a function of energy (Fig. [2]) . The average event rate 
of electron recoil singles below 100 keV in all detectors is 
stable in time at 16%. 

For the germanium detectors considered in this anal- 
ysis, the summed background rate after correcting for 
detection efficiency is ~1.5 cpd (counts per day) kg" 1 
keV" 1 (Fig. [3]). The prominent 10.36 keV line is caused 
by X-rays and Auger-electrons from the electron-capture 
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FIG. 2: Detection efficiency as a function of energy. 



decay of 71 Ge, produced by neutron capture on 70 Ge 
during 252 Cf calibration of the detectors. The excess in 
event rate around 6.5 keV (inset) is likely caused by rem- 
nant 55 Fe decays from cosmogenic activation. The de- 
excitation of 55 Mn following the electron-capture decay 
of 55 Fe yields a total of 6.54 keV of electron-recoil events. 
We interpolate the energy resolution of the 10.36 keV line 
(typically better than 3%) to the noise level to obtain the 
energy-dependent resolution of each detector. The anal- 
ysis window, defined from 2 to 8.5 keV, is determined by 
the expected axion flux, background rate, and detection 
efficiency. 

We carried out a profile likelihood analysis to deter- 
mine the best fit value of 5077- We express the event rate 
per unit measured energy (E), per unit time (t), and per 
detector (g0 of a solar axion signal with background as 



R(E, t, d) = s{E, d) [\K(E, t, d) + B(E, d)} 



(3) 



where e(E,d) is the detection efficiency, lZ(E,t,d) is 
the expected event rate for a coupling constant 5,377 = 
lfT 8 GeV -1 , and A = (g ail -10 s GeV) 4 is the scale factor 
for the actual value of g a ~n- B{E 1 d) is the background 
described by 



B{E,d) = C{d) + D(d)E + H{d)/E 



Vs. 54 



(E-6.54kcV) i 





(4) 



54 



where C(d), D{d), and H{d) are free parameters. The 
Gaussian term describes a contribution from 55 Fe decays 
at an energy of 6.54 keV and unknown total rate 776.54- 
The fitting is done by maximizing the unbinned log like- 
lihood function with respect to A and the background 
parameters, for individual events i: 

log(C) = -R T + V log(R(Ei, U,dj)), (5) 

where Rt is the total sum of the event rate (R) over en- 
ergy, time, and detectors. The scaling factor from the 
maximization A = (1 ± 1.5) x 10~ 3 is compatible with 
zero. No indication of solar axion conversion to photons 
is observed. Given a null observation, we set an upper 
limit on the coupling constant ga-yj, where the scaling fac- 
tor A is obtained by integrating the profile likelihood in 




FIG. 3: Co-added, efficiency corrected low energy spectrum of 
the Ge detectors considered in this analysis. The inset shows 
an enlargement of the spectrum in the analysis window, taken 
to be 2 to 8.5 keV. 



the physically allowed region (A > 0). The upper limit on 
the axion-photon coupling, g ail < 2.4 x 10~ 9 GeV -1 at a 
95% C.L. is the only laboratory bound based on the accu- 
rate measurement of all crystal orientations of the detec- 
tors. None of the previous solar axion search experiments 
(SOLAX/COSME/DAMA) measured their crystal orien- 
tations I15l . , and thus their limits are penalized by 
picking the least sensitive direction for their solar axion 
bound. The result of this analysis is compared to other 
experimental constraints in Fig. [5] Improvement towards 
the next order of sensitivity requires improvements in 
both detector exposure and gamma background level. A 
100-kg SuperCDMS experiment, with substantially re- 
duced gamma background level (~0.1 cpd kg -1 keV -1 ) 
would improve the sensitivity to g ail < 10 -9 GeV -1 . 

In addition to restricting solar axions, the CDMS mea- 
surement can be used to limit Galactic axions. The 
annual modulation signature observed by DAMA may 
be interpreted as a detection of axions distributed in 
the local Galactic halo [2^, [2l[. If present, these non- 
rclativistic axions would materialize in our detectors via 
an axio-electric coupling (g a se)- Assuming a local Galac- 
tic dark matter mass density of 0.3 GeV/c 2 /cm 3 , the ex- 
pected event rate [22| is given by: 



R [cpd kg -1 ] = 1.2 x 10 43 A- 1 g 2 aSe m a a p 



(6) 



where m a is the axion mass in keV/c 2 , ^4=73 for ger- 
manium, and (jp.e is the photoelectric cross section in 
cm 2 per atom. We analyzed the energy spectrum using 
the same electron-recoil data samples used in the solar 
axion search, as shown in Fig. [3l We carried out a pro- 
file likelihood calculation to search for an excess of event 
rate above background. The same formalism described 
in equations ([3]) to §5§ was used, with the term for the ex- 
pected conversion rate of solar axions lZ(E,t,d) replaced 
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FIG. 4: Comparison of the 95% C.L. upper limit on 
g a ~/~i achieved in this analysis (red/solid) with other 
crystal search experiments (SOLAX/COSME [H, Q]| 
(black/solid) and DAMA (upper black /dashe d) PulO and he- 
lioscopes (Tokyo heslioscope (magenta/solid) jl8| ] and CAST 
(blue/solid) [llj]). The constraint from Horizontal Branch 
stars (lower black/dashed) is also shown [ljj. 
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FIG. 5: The allowed region (green/filled) from a Galactic 
axion interpretation of the annual modulation signature ob- 
served by the DAMA experiment is shown. The 90% C.L. 
upper limits on the g a ee coupling constant from this work 
(red/solid) and the CoGeNT experiment (blue/solid) [23j] 
completely exclude the DAMA allowed region. The indirect 
constraints from astrophysical bounds (black/dashed) are also 
shown [24[. 



by a Gaussian distribution function representing a spec- 
tral line at a given energy or axion mass. We find no 
statistically significant excess of event rate above back- 
ground. Lacking a direct constraint on a possible 55 Fe 
contribution to the spectrum, we set a conservative up- 
per limit, shown in Fig. [5j on the total counting rate in 
this energy range without any attempt to subtract a pos- 
sible background contribution. This result excludes sig- 
nificant new Galactic axion parameter space in the mass 
range between 1.4 and 9keV/c 2 , and is inconsistent with 
the interpretation of the DAMA signature due to axions 

In summary, the solar axion search sets an upper limit 
on the Primakov coupling g ail of 2.4x 1CP 9 GeV -1 at 
the 95% confidence level for an axion mass less than 
^O.lkeV/c 2 . This limit is the first one based on accu- 
rate measurements of crystal orientations. The system- 
atic error on the limit is estimated to be 7.9%, which 
arises from the remaining uncertainty in the alignment 
of the detector towers' major axes to the central cryo- 
genic axis. The local Galactic axion search analysis sets 
a world-leading experimental upper limit on the axio- 
electric coupling g aSe of 1.4x 10 -12 at the 90% confidence 
level for an axion mass of 2.5keV/c 2 , and excludes the 
axion interpretation of the DAMA signature in the mass 
range of 1.4keV/c 2 to 9keV/c 2 . 

This experiment would not have been possible without 
the contributions of numerous engineers and technicians; 
we would like to especially thank Larry Novak, Richard 
Schmitt and Astrid Tomada. We thank the CAST and 
Tokyo helioscope collaborations for providing us with 
their axion limits. The direction measurement of the true 



north in the Soudan Underground Laboratory relied on 
the help from the Fermilab Alignment Group. Special 
thanks to Virgil Bocean. This work is supported in part 
by the National Science Foundation (Grant Nos. AST- 
9978911, PHY-0542066, PHY-0503729, PHY-0503629, 
PHY-0503641, PHY-0504224 and PHY-0705052), by the 
Department of Energy (Contracts DE-AC03-76SF00098, 
DE-FG02-91ER40688, DE-FG02-92ER40701, DE-FG03- 
90ER40569, and DE-FG03-91ER40618), by the Swiss 
National Foundation (SNF Grant No. 20-118119), and 
by NSERC Canada (Grant SAPIN 341314-07). 



[1] R. D. Peccei & H. R. Quinn, Phys. Rev. Lett. 38, 1440 
(1977). 

[2] S. Weinberg, Phys. Rev. Lett. 40, 223 (1978). 

[3] F. Wilczek, Phys. Rev. Lett. 40, 279 (1978). 

[4] J. E. Kim, Phys. Rev. Lett. 43, 103 (1979); J. E. Kim, 

Phys. Rep. 150, 1 (1987). 
[5] P. Sikivie, Phys. Rev. Lett. 51, 1415 (1983). 
[6] K. D. Irwin et al, Rev. Sci. Instr. 66, 5322 (1995); 

T. Saab et al, AIP Proc. 605, 497 (2002). 
[7] D. S. Akerib et al, Phys. Rev. D 72, 052009 (2005). 
[8] Z. Ahmed et al, Phys. Rev. Lett. 102, 011301 (2009). 
[9] J. Bahcall & M. H. Pinsonneault, Phys. Rev. Lett. 92, 

121301 (2004). 

[10] S. Andriamonje et al, J. Cosmol. Astropart. Phys. 04, 

010 (2007). 

[11] E. Arik et al. larXiv:0810.4482l 

[12] E. A. Pascos & K. Zioutas, Phys. Lett. B 323, 367 (1994). 

[13] S. Cebrian et al, Astropart. Phys. 10, 397 (1999). 

[14] Fermilab Alignment Group survey at Soudan Under- 



5 



ground Laboratory, drawing H5.1, MN1002/3 (1999). 
[15] F. T. Avigone III et ai, Phys. Rev. Lett. 81, 5068 (1998). 
[16] A. Morales et al, Astropart. Phys. 16, 325 (2002). 
[17] R. Bernabei et ai, Phys. Lett. B 515, 6 (2001). 
[18] M. Minowa et al, Phys. Lett. B 668, 93 (2008). 
[19] G. G. Raffelt, Stars as Laboratories for Fundamental 

Physics, The University of Chicago Press (1996). 
[20] R. Bernabei et ai, Int. J. Mod. Phys. A 21, 1445 (2006). 
[21] R. Bernabei et al, Eur. Phys. J. C 56, 333 (2008). 
[22] M. Pospelov, A. Ritz & M. Voloshin, Phys. Rev. D 78, 

115012 (2008). 



[23] C. E. Aalseth et al., Phys. R ev. Lett . 101, 251301 (2008). 
[24] P. Gondolo & G. G. Raffelt. larXiv:0807.2926l 
[25] The DAMA region should be understood with caution. 
The non-relativistic speed of Galactic axions distributed 
in the local halo causes the conversion rate to be indepen- 
dent of the particle's velocity, thus the annual modulation 
of the counting rate is highly suppressed [22|. Therefore, 
the DAMA modulation signal is too large to be inter- 
preted by Galactic axion interactions in the detectors. 



